Interference effects in the reradiation of ultrashort electromagnetic pulses 
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Reradiation of a spatially non-uniform ultrashort electromagnetic pulse interacting with the linear 
chain of multielectron atoms is studied in the framework of sudden perturbation approximation. 
Angular distributions of the reradiation spectrum for arbitrary number of atoms are obtained. 
It is shown that interference effects for the photon radiation amplitudes lead to appearing of 
"diffraction "maximums. The obtained results can be extended to the case of two- and three- 
dimensional crystal lattices and atomic chains. The approach developed allows also to take into 
account thermal vibrations of the lattice atoms. 
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Interaction of short and ultrashort electromagnetic 
pulses with atomic systems and discrete structures 
such as crystal lattices, atomic and molecular chains 
is of importance both from practical and fundamental 
viewpoints. Such processes can be used in X-ray structure 
analysis and attosecond spectroscopy increasing their 
resolution. Attosecond physics uses as the targets 
relatively simple systems and their similar extensions Q- 
[5] in the context of ultrashort processes During 
recent years one can observe rapidly growing interest 
to the interactions of ultrashort pulses which is caused 
by their potential applications in supestrong lasers and 
heavy ions accelerators where strong fields of highly 
charged ultrarelatuvistic ions are comparable with those 
of light field. In particular, in the collision experiments of 
with fast highly ions the field of ultrarelativistic uranium 
ion, U 92+ with the collision energy 1 GeV/nucleon can 
be considered as superintense pulse with the strength 
(I > 10 19 Wt/cm 2 )and the duration r ~ 1CT 18 s. 
Despite certain progress made in the study of interactions 
of ultrashort pulses with atomic systems, reradiation 
induced by such interaction is still remaining as less- 
explored topic which is considered only few papers (see. 
e.g. Q and references theirein). 

Usually diffraction of X-rays on different periodic 
structures is described as the plane wave scattering 
with infinite time duration pj. In Q the scattering 
of ultrashort pulse on the atom is studied within 
the classical approach. The scattering of ultrashort 
electromagnetic pulse in multielectron atoms within the 
mantum mechanical approach is considered in the Refs. 



quant 
0,1. 



In particular, in 0, Q the spectra of scattered 
radiation are obtained for different values of the 
pulse duration. Such an approach can be applied for 
attosecond pulses, too. However, exact treatment for 



such ultrashort pulses is possible within the sudden 
perturbation approach which enables also the treatment 
of the re-scattering processes [IcJ and can be extended 
to the cases of simplest molecules 

Below we consider ultrashort pulses whose duration, r 
is much shorter that the characteristic atomic period r a 
for the Bohr orbit. 

In th is p aper using the approach developed in the 
Ref. [10(, we treat reradiation of the ultrashort 
pulses of the electromagnetic field in the liner chain 
of the multielectron atoms. The method developed 
allows exact account of the non-uniform distribution 
of both ultrashort pulse as well as photon momentum 
for reradiation processes. The field of the ultrashort 
pulse can be taken into account within the sudden 
perturbation, while photon radiation is described by the 
usual perturbation theory. Also we calculated angular 
distributions of the reradiation spectra for arbitrary 
number of atoms in the chain. It is shown that 
interference effects in the amplitudes of the radiated 
photons may lead to appearing of "diffraction "maximums 
which become infinitely narrow in the limit of infinitely 
large number of atoms. 

Consider an atom interacting with the electromagnetic 
pulse of Gaussian shape, (atomic units are used 
throughout the paper): 

E(r, t) = E ^e- a2( *- nor/c)2 cob(w * - k r) , (1) 

\/7T 



where k = (a; /c)n , a n is the unit vector directed 
along the pulse direction , r is the coordinate of the 
observation point, c is the light speed and r ~ 1/ot is the 
pulse duration. We note that E(r, t) — > E <5(£ — n r/c) 
for a — > oo. 

According to the Ref. [10j the interaction of the atomic 
electrons with electromagnetic pulse can be described by 
the potential 
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V(t) = V({r e },t)= E(r e ,i)r e 



(2) 
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with {r e } being the collective coordinate of the atomic 
electrons and (e = 1, N e ), N e is their number. 

Let the value of the quantity a in Eq. ( [T]) such that 
V(t) is effectively non-zero only during the time, r ~ 
a -1 , which is much shorter that the characteristic periods 
of atomic electrons whose dynamics is described by the 
Hamiltonian Hq. 

Then the amplitude for transition of the atom from the 
initial , ip top a final state, tp n caused by the interaction 
with sudden perturbation, V(t) can be written as (l2| : 



+00 



aon = (Pn I exp(~i / V(t)dt) | <p ), 



(3) 



where ip a and ip n are the eigenfunctions of the 
unperturbed Hamiltonian, H . Perturbation potential 
given by Eq. ( [5]) is written for the atom whose nucleus is 
at the origin of coordinate. For the atom with the nucleus 
at the distance d from the origin of coordinate system the 
interaction potential can be written as ( [5]) 

e—N e 

V(t) = E(r e +d,t)(r e +d) = 

e=l 

e=IV e 

= E ( r e + d > *) r e + E(r e + d, t)d , (4) 



e=l 



where r e are the distance between the atomic electrons 
and nucleus. It follows from Eq.([T]) that for E(r e + d,t) 

+00 

the integral J dt from second term in the right hand side 

—00 

of Eq.(0| ) doesn't depend on the coordinates of atomic 
electrons. Therefore it doesn't make contribution into the 
transition amplitude, ao n . 

Thus within the sudden perturbation approach as the 
interaction potential one can use 



V(t) = V(r e ,t) = E(r e + d,i)r e 



(5) 



Now consider one-dimensional chain of N 
noninteracting atoms with the same distances, d 
between neighboring atoms and denote by N e the 
number of electrons in each atom. Furthermore, we 
assume that the origin of coordinate coincides with 
the position of the first atom, so that the distance 
between the first and a— th atoms is (a — l)d. Then 
the coordinate of the electron in a— th atom is given 
by R a!e = (a — l)d + r ayC , where r 0jE is the distance 
between the electron and nucleus in the a— th atom 
(o = l,2,---,JV). 

Using these notation the potential for interaction of the 
chain with a pulse of an external electromagnetic field can 
be written as 



N N e 



v{t) = J2Y1 E ( r <^ + ( a - 1 ) d ' 



(6) 



To consider this potential as a sudden perturbation for 
a single atom the following condition should be obeyed: 
r ~ 1/a -C r a . However, for a chain of atoms with the 
length L — N\d\ the potential V(t) can be considered as 
a sudden perturbation under the following conditions: 



r — 1/a < T a ~ 1, 
iVdn /c< r a - f. 



(7) 



We assume that each atom in the chain initially being 
the ground state, ip ({r aye }), with {r a!e } being the set of 
coordinates of the atomic electrons. Then for the wave 
function of the one-dimensional chain consisting of N 
noninteracting atoms with total NN e electrons we have 



$0 = <A)({ri, e })<A)({r2, e }) • • • <p ({r Nie }). 



(8) 



In the following, we will be interested in the treatment 
of reradiation of the ultrashort electromagnetic pulse 
during the time its interaction with the above atomic 
chain. 

Following the Ref. [10| the interaction of the atomic 
electrons with the electromagnetic field can be written as 



N N e 



a—1 e—1 k,cr 



+ -ikR a cA 

Uk CT a kcr e ac p a 



(9) 



Here a^ a are the photon creation operators, u), k and 
a, (a = 1,2) are the frequency, momentum and 
polarization of the created photon, respectively, Uko- are 
the unit vectors of polarization, while p a e = d/dr a ^ e are 
the momentum operators of the atomic electrons. Then 
the amplitude for the photon emission can be calculated 
in the first order of perturbation theory with respect to 
interaction U [loj . 

We want to find the photon radiation spectrum of on 
the solid angle dfl^ , drawn by the vector k. For this 
purpose we use an approach developed in the Ref. flfj| 
which gives us (averaged over photon polarization vector 
and summed over all the final states of atoms in the 
chain) the total spectrum of the photon radiation per 
unit of solid angle dfl^: 



d 2 W 



dQ^duj (27r) 2 c 3 lj 



a,e b,m 



-ik(Ro, e — R(,,„ 



xf^-^nH 7 n] I $ ),(10) 



dr a . 



dr b . 



where R Qie = (a— l)d+r a)e and Rh e = (b— l)d + r(,. e are 
the coordinates of the electrons belonging to a— th 6— th 
atoms, respectively, n = k/fc is the unit vector along 
the direction of emitted photon and V(u>) is the Fourier- 
transform of the potential V(t), which is written as 

+00 

V(lj)= J V{t)e iut dt = f (Lj)Y,(Vora,e)x 



xexpli — k ((a- l)d + r . e ) , (11) 
w J 
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with 



/oM = -< exp 



+exp 



(u) - uj ) 




(12) 



and 



dVioS) { uj 1 

[— n] = f (uj)exp < i — ko((a - l)d + r a , e ) ^ x 

x (jEon] + i-^(Eor a , e )[konf) .(13) 

Using Eqs. (ITTj) - p^l the total radiation spectrum per 
unit od solid angle dQ^ for one pulse duration, r can be 
written as 



d 2 W | / (w) | 5 



dfik^ (27r) 2 c 3 w 



6) 



-ip(r...-r k , m ) J [ Eo n] 2 +i — ([E n][k n]) x 
xEo(r a!e - r b . m ) + 



-(E r , e )(Eor 6 . m )- ? [k n] 2 } | $„). (14) 



It is clear that in this formula, the averages over 
the ground state, <&o are written in the form (<E> | 
X) f( r a.e,Yb,m) \ &o), and therefore do not depend on 

e,m 

a and b. 

Separating in Eq. (|14[) summation for a — b and a ^= b 
and denoting corresponding spectra by c? 2 Wi/(o?ilkrfw) 
and d 2 W2/(dfikrfw), respectively, we have 



d 2 W d 2 W x d 2 W 2 



dflkduj dQkduj dil^duj 



(15) 



In the following we will see that the term 
d 2 Wi j (dClkdu) describes incoherent (proportional to 
N) rescattering , while the interference spectrum 
d W2/(dQ.kdu), corresponds to appearing of 
characteristic diffraction maximums. 

Let us apply the above formalism to the linear atomic 
chain consisting of N non-interacting helium atoms. If 
yo( r ii r 2) is the ground state wave function of helium 
atom, then the wave function of the chain, $o which 
appears in Eq. (|14[) can be written as 

$o = ¥ ? Q( r i,i ! r i,2)^o(r2,i ! r 2 ,2) ■ • • <Po(tn,i,*n,2)- 

Furthermore, we write each ground state wave function 
of the helium atom, (p as the product of two hydrogen 
like Is wave functions with an effective charge Z = 2 — 



5/16. Then one can easily perform averaging over $o hi 
Eq. (Ti"4")) and the spectra can be written as 

+N e (N e -l)F{Lj,n,n )\, (16) 



d 2 W 2 | /„(«) | 5 



dQkdiu (2tt) 2 c 3 uj 



7V 2 F( W ,n,n ) ffJV (pd) , (17) 



where 



G(«,n,no) = [E n] 2 



F(u,n,n ) 



c 2 Z 2 
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£; 2 [n n] 2 , (18) 



[4+(n-n ) 2 o; 2 /(cZ) 2 ] 2 



x <^ [E n] 



n n 



4(E (n-no)) 



c 2 Z 21 u J 4+(n-n ) 2 w 2 /(cZ) 2 



g N {pd) = e 

pd 



pd(a-b) _ sin 2 (pdA^/2) _ 



sin 2 (pd/2) 



,(19) 



N, (20) 



(n - n )d. 



(21) 



Here n and no are the unit vectors along the directions 
of the emitted photon and incident ultrashort pulse, 
respectively, oj is the frequency of the emitted photon, 
d is the vector directed along the chain such that |d| is 
the distance between the neigbouring atoms. Summing 
Eqs. (HHJ) and (jTTJ) we get the total radiation spectrum 



d 2 W 



dflkduj (2it) 2 c 3 uj 
-F(u,n, no) 



G(w,n, n ) 



„.„, N sin 2 (pdA^/2)l 
7V e 2 F(^,n,n ) , f > \. 

sm (pd/2) 



(22) 



This equation is valid for a chain with arbitrary N 
provided the condition given by Eq. ([7])is obeyed. For 
Z = 1 h N e = 1 Eq. (|22p describes reradiation of the 
chain of N hydrogen atoms induced by the interaction 
of the chain with an ultrashort electromagnetic pulse. 
For N — 1 and N e = 1 it describes reradiation of single 
hydrogen atom while N = 1 and N e = 2 corresponds to 
the reradiation of single helium atom (see (Io|). 

It is easy to see that near the zeros of its denominator, 
second term in right hand side of Eq. (|22|) behaves like 
iV 2 , while for for N ^> 1 it reproduces the behavior 
which is usual for of the angular distribution of the wave 
processes on diffraction gratings. It should be noticed 



4 



that the functions G(ui, n, no) h F(ui, n, no), contain (see 
Eqs. (fT5|) and (|19p L the frequency of radiated photon, u> 
only in the form uj 2 /(cZ) 2 . 

Now let us consider behavior of of the spectrum given 
by Eq. (|2"2")l in the limit of small frequencies given 
by inequality ui 2 /(cZ) 2 <C 1. In this limit we have 
G(w, n,n ) = F(w,n, n ) = [E n] 2 and Eq. can 
be written as 



d 2 W 



(irik^w (27r) 2 c 3 w 



7V e 2 [E n] 



2 sin 2 (pdA^/2) 
sin 2 (pd/2) 



(23) 



The contribution from the interference effects can be 
characterized by the quantity 



d-*W 

d 2 W t 
dQ\rdLj 



1 sin 2 (pdiV/2) 
sin 2 (pd/2) 



(24) 



In this formula zeros of denominator in the right hand 
side correspond to maximum values of the angular 
distribution. In the limit N ^> 1 for pd/2 = irn + e (n is 
integer and, e is a small quantity) , we have an asymptotic 
relation [13j 

hm l**Wm = Um 1^50=^25) 
n^oo N sin 2 (pd/2) JV-s-oo N e 2 

It is easy to see that in Eq. (|22[) near the above 
mentioned zeros of denominator second term behaves 
like (N e N) 2 and corresponds to coherent radiation of 



photon by all electrons in the chain. As it can been from 
Eq.^for u 2 /{cZ) 2 ->■ oo (or, equivalently, uj 2 /{cZ) 2 > 
NN e ) the coherent part of the spectrum is much smaller 
compared to incoherent part (which is proportional to 
NN e ). Therefore increasing of to leads to decreasing of 
the contribution by coherent radiation. 



Eq. (|22p can be utilized for qualitative description 
of reradiation process caused by the interaction with 
spatially non-uniform ultrashort pulse for arbitrary 
number of atomic electrons (N e > 2) and atoms in 
the chain. It follows from Eq. (fl"4"]) that the dependence 
of the total spectrum on N and N e has the same 
character as that in Eq. (|22D . The only difference comes 
from the functions G(uj, n, no) and F(u>, n, no) which 
are calculated as the averages over the ground state 
wave functions of some atoms in the chain. In the 
above treatment we considered atoms in the chain as 
noninteracting to each other ignoring by electron-electron 
correlations. However, screening of nucleus charge by 
inner electrons has been taken into account using the 
effective charge in Eqs. (TT51) and (TT71) . Finally we 
note that the above treatment of reradiation process 
induced in linear atomic chain by its interaction with 
a spatially non-uniform ultrashort pulse can be extended 
to the cases of higher dimensional (e.g., graphene, three- 
dimensional crystall lattices etc) atomic chains, too. 
Also, it is possible to include into consideration thermal 
vibrations of the atoms in the chain. Such a study is 
currently on progress. 
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